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INTRODUCTION 


This  report  discusses  the  use  of  adjoint  variable  formulation  to  seek  the 
transient  solutions  for  problems  in  gun  dynamics.  The  theory  from  variational 
principle  involving  adjoint  variables  solves  a  mixed  boundary  and  initial 
value  problem.  The  partial  differential  equation  governing  the  motion  has  a 
fourth  order  partial  in  spatial  domain  and  a  second  order  partial  in  time 
domain.  It  also  involves  a  few  step  functions  and  delta  functions  as 
follows . 1 >  2 

pAy  +  (Ely”)*’  -  [PCx^y*]'  +  Ty"H(x-xp)  - 
ra[xp2y"  +  2xpy '  +  y]6(x-xp) 

-  mg  cos  a  S(x-Xp)  -  pAg  cos  a  (1) 

The  above  equation  can  be  simplified  into  the  following  form 

Ly  +  Q  *  0  (2) 

where 

Ly  =  (ayt)t  -  (^Yxx)xx  +  (^Yx)x  +  (^pYx)x^(x“xp)  (3) 

and 

-Q  =  m[xpZy"  +  2xpyf  +  y]S(x-xp)  -  mg  cos  a  6(x-xp)  -  pAg  cos  a  (4) 

We  seek  the  explicit  numerical  transient  solutions  of  y,  yt ,  yx,  Yxt >  Yxx > 
and  yXxt  ^or  sorae  given  boundary  arid  initial  conditions.  The  terra  yxx  will 
give  the  stress  wave  and  the  terra  yx  will  show  the  slope  in  bending,  along 


ISimkins,  Thomas  E. ,  "Transverse  Response  of  Gun  Tubes  to  Curvature-Induced 
Load  Functions,"  presented  at  the  Second  US  Army  Symposium  on  Gun  Dynamics, 
Watervliet,  NY,  September  1978. 

^Wu,  Julian,  "Gun  Dynamic  Analysis  by  the  Use  of  Unconstrained,  Adjoint 
Variational  Formulations,”  presented  at  the  Second  US  Army  Symposium  on  Gun 
Dynamics,  Watervliet,  NY,  September  1978. 
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the  axis  of  the  gun  tube.  The  solution  is  the  extension  of  our  previous 
work  on  initial  and  boundary  problems. 

VARIATIONAL  PRINCIPLE  USING  ADJOINT  VARIABLE 

If  the  inner  product  of  the  variable  y,  and  the  adjoint  forcing  function 
Q  are  used  for  variational  purposes,  the  accuracy  is  much  less  than  the  method 
using  the  following  inner  product  by  adding  a  terra  involving  the  adjoint 
variable  y  as  the  Lagrange  multiplier  (see  Appendix). 

Jty.y]  =  <Q»y>  +  <y,(Q+-Ly)>  =  0  (5) 

where  the  partial  differential  equation  is  given  in  Eq.  (2).  By  taking 
variation  on  Eq.  (5)  we  have 

SJ  =  <Sy,(Ly+Q)>  +  <Sy,(Ly+Q)>  -  <6y,Ly>  +  <y,L6y>  =  0  (6) 

The  above  variation  vanishes  if 

Ly  +  Q  =  0  (7) 

Ly  +  Q  =  0  (8) 

and  _  _ 

<y»L6y>  -  <6y,Ly>  =  0  (9) 

We  know  that  Eq.  (7)  is  actually  the  original  p.d.e.  and  Eq.  (8)  is  its 
adjoint  equation.  A  method  should  be  established  so  that  Eq.  (9)  holds  true 
for  all  arbitrary  variation  6y. 


JShen,  C.  N.  and  Wu,  J.  J. ,  "A  New  Variational  Method  for  Initial  Value 
Problems  Using  Piecewise  Hermite  Polynomial  Spline  Functions,"  presented  at 
the  1981  Army  Numerical  Analysis  and  Computer  Conference,  Huntsville,  AL, 
February  1981. 

4Shen,  C.  N. ,  "Method  of  Solution  for  Variational  Principle  Using  Bicubic 
Hermite  Polynomial,"  presented  at  the  27th  Conference  of  Army  Mathematicians, 
West  Point,  NY,  June  1981. 
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BILINE  Ail  CONCOMITANT 


We  will  find  out  the  conditions  for  the  assumed  equality  in  Eq.  (9)  to  be 
true.  Let  us  consider  the  following  bilinear  concomitant : 5 

D  -  <y,Ly>  -  <y,Ly>  (10) 

The  above  expression  can  be  integrated  in  two  different  ways  and  can  also 
be  written  in  terms  of  boundary  conditions  and  initial  conditions.  It  is 
assumed  that  these  boundary  conditions  are  assigned  in  such  a  manner  that  the 
above  bilinear  concomitant  is  identically  zero  for  all  independent  variables, 
i.e.  , 


o 

Ml 

Q 

(11) 

Then  the  fir9t  variations  of  D  also  vanish. 

6D  =  6D(6y)  +  6D(6y)  =  0 

(12) 

Since  6y  and  6y  are  independent  of  each  other,  then 

6D(6y)  =  <6y ,Ly>  -  <y,L6y>  =  0 

(13) 

SD( 6y)  =  <y,L6y>  -  <6y,Ly>  =  0 

(14) 

Equation  (14)  is  identical  to  Eq .  (9),  which  is  the  assumed  equality 
previously.  The  implication  is  that  if  Eq.  (11)  is  true  then  Eq.  (9)  or  (14) 
is  automatically  true. 

Since  Eq.  (10)  can  be  expressed  in  terms  of  some  integrals  involving 
boundary  conditions,  Eq.  (11)  can  be  true  if  these  boundary  conditions  are 
satisfied.  The  next  section  will  discuss  integral  of  bilinear  expression  and 
its  boundary  conditions. 

^Stacey,  W.  M.  Jr.,  Variational  Methods  in  Nuclear  Reactor  Physics ,  Academic 
Press,  1974. 
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INTEGRAL  OF  BILINEAR  EXPRESSION 


The  integral  of  a  bilinear  expression  for  a  two  dimensional  problem 

having  second  order  partial  derivatives  in  time  and  fourth  order  partial 

derivatives  in  space  can  be  written  as 

xb  tb 

1  =  /  /  U[y(x,t)y(x,t)]dtdx  (15) 

xo  to 

where  ft[y,y]  is  a  given  bilinear  expression  in  the  form 

nly»y]  =  «ytyt  +  Xyxxyxx  +  *yxyx  +  **yxyxH(x-xp)  06) 

The  subscripts  t  and  x  indicate  the  partial  derivatives  of  the  functions  y  and 


Equation  (16)  can  be  integrated  by  parts.  Two  different  forms  of 

integration  and  end  conditions  are  given.  The  first  form  of  the  integral  is 

obtained  by  integrating  by  parts  on  the  adjoint  variable. 

tb  xb  -  xb  tb 

I  ■  -/  /  yLydtdx  +  /  cxyty  |  dx  + 

to  xo  xo  t0 

tb  -  xb  -  xb  xp  xb 

I  t^yxxYx I  “  (^yxx)xy|  +  ^yxy|  +  ^*yxy|  Jdt  (17) 

to  xo  xo  xo  xp 

where 

Ly  =  (ayt)t  ~  Oyxx)xx  +  Oyx)x  +  0*yx)xH(x-xp)  08) 

On  the  other  hand,  we  can  perform  integration  on  the  original  variable  to 

give 

tb  xb  __  xb  -  tb 

I  =  “J  j  yLydtdX  +  /  ayty|  dx 
t0  xo  xo  tQ 

tb  _  xb  -  xb  _  xb  Xb 

+  J  Oyxxyx |  -  Oyxx)xy |  +  £yxy|  +  ^*yxy|  }dt  (19) 

to  xO  xO  xO  xp 
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where 


Ly  =  (ayt)t  -  (Xyxx)xx  f  Uyx)x  +  ( £pyx)xH(x-xa)  (20) 

For  a  fourth  order  spatial  partial  and  a  .second  order  temporal  partial  system 
Eq.  (10)  becomes 

xb  tb  -  xb  tb  — 

D  =  /  J  yLydtdx  -  J  J  yLydtdx  (21) 

xO  to  X0  to 

By  equating  Eqs.  (17)  and  (19)  and  solving  for  D  in  Eq.  (21)  we  are  converting 
the  double  integral  into  two  single  integrals  in  terms  of  the  boundary 
conditions. 

We  can  express  the  quantity  D  as  the  sum  of  three  parts  on  end  conditions 
D^,  D2,  and  D3  as 


D  =  Dj  +  D2  +  D3  '  (22) 

The  terms  in  involve  the  initial  conditions  of  y  and  y  as 

xb  -  tb  tb 

Di  =  /  Wtyl  -  “ytyl  ^dx 

xo  tQ  tQ 

xb  -  -  -  - 

=  J  (“b(ytbyb-ytbyb)  -  «o(ytoyo~ytoyo)>dx  (23) 

xo 

The  terms  in  D2  involve  the  boundary  conditions  from  the  second  partials  of  y 
and  y  as 


tb  xb  -  xb  -  xb  -  xb 

D2  =  /  Uyxy|  -  *yxy|  +  ^yxy  I  -  **yxy|  >dt 

to  xo  XQ  Xp  Xp 


-  /  Ubyxbyb  ~  ^oyXoyo  -  ^byxbyb  +  ^oyxoyo 
to 

+  ^b*yxbyb  -  £PyxPyP  -  *byXbyb  +  Vxpyp^t  (24) 
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The  terms  in  D3  involve  the  boundary  conditions  from  the  fourth  partials  of  y 
and  y  as 

fcb  -  *b  -  *b  -  xb  _  xb 

“  J  Uyxxyx I  —  (^yxx)xy|  -  ^yxxyxl  +  (^yxx)xyl  Jdt 

C0  *0  Xq  Xq  Xq 

fcb 

=  J  ^byxxbyxb  -  *oyxxoyXo  -  (^yxx)xbyb  **■  (^yxx)xoyo 
co 

-  ^byxxbyxb  +  ^oyxxoyxo  +  (^yxx)xbyb  -  (^yxx)xoyo^t  (25) 

In  order  that  D  =  0  in  Eq.  (22)  it  is  sufficient  that 

Di  =  0  (26a) 

D2  =  0  (26b) 

and  D3  =  0  (26c) 


END  CONDITIONS  FOR  THE  ADJOINT  SYSTEMS 

In  order  to  satisfy  the  requirements  in  Eq.  (26)  we  separate  them  again 
in  three  different  parts. 

(a)  Let  us  assume  that  the  adjoint  variables  are 

yb  =  kiyo  »  y0  =  kiyb  (27) 

ytb  =  -ab_laokiyto  »  yto  =  -ao_labkiytb  (28) 

where  kj  is  a  constant.  The  above  adjoint  boundary  conditions  satisfy  the 
requirement  that  Dj  =  0  in  Eq.  (23). 

(b)  Let  us  assume  the  following  adjoint  variables 


yb  =  k2yb 


yxb  ”  k2yxb 


-  kl2 

y0  =  , —  y0 

k2 


kl2 

yxo  =  — -  yxo 

k2 


yp  -  k3yp 

yxp  =  k3yxp 


(29) 


(30) 
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Where  Eq.  (29)  Is  Inconsistent  with  Eq .  (27)  and  k2  is  another  constant. 
Equations  (29)  and  (30)  Imply  that  D2  =  0  in  Eq.  (24). 

(c)  The  following  boundary  conditions  are  assumed 

_  k^2  _  _  -  k}2 

To  =  Yo  j  Yxo  “  Yxo  >  Yxxo  =  Yxxo  >  Yxxxo  “  (“  )  Yxxxo  (31) 

k2  k2 

Yb  =  ^2Yb  >  Yxb  =  Yxb  >  Yxxb  =  Yxxb  >  Yxxxb  =  k-2  Yxxxb  (32) 
Equations  (31)  and  (32)  satisfy  Eq.  (25).  Thus  D3  =  0. 

By  giving  the  appropriate  values  of  the  adjoint  variables  in  terms  of  the 
original  variables  one  may  find  that  the  requirement  D  =  0  can  be  satisfied. 
This  leads  to  the  condition  in  Eq.  (6)  that 

5J  =  0 

for  all  arbitrary  variations  5y  and  5y. 


FIRST  VARIATION 

Since  the  variations  6y  and  6y  are  independent  to  each  other,  the  part  of 
6j  in  Eq.  (6)  with  variation  6y  can  be  expressed  as 

xb  tb  -  xb  tb  - 

6j(6y)  =  /  /  5yLydt  +  /  /  6yQdtdx  =  0  (33) 

xo  to  xo  tQ 

Where  Ly  is  given  in  Eq.  (18)  and  contains  second  and  fourth  partial 
differentiations  in  y.  It  is  intended  to  include  only  low  order  partial 
differentiations  in  6J(5y).  This  can  be  achieved  by  considering  the 
variations  of  the  bilinear  expression  I  given  by  Eqs.  (15)  and  (16)  as, 

tb  xb 

6j(6y)  =  J  J  [ayt<5yt  +  ^xx^xx  +  S-yxSyxHtdx 


tb  xb 

+  J  J  ^pyx'5yx<it;dx 
to  x0 


(34) 
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A  different  form  of  the  above  variation  can  be  obtained  from  Eq.  (17)  as 


xb  -  tb 

6l(6y)  =  -//  6yLydtdx  +  /  ayt6y|  dx 

xo  ^o 

tb  -  xb  -  Xb  -  Xb  -  xb 

+  /  ^yxx^yxl  -  (*yxx)x6y|  +  *yx5y|  +  £*yx6y|  ^dt  (35) 

tQ  X0  Xo  xo  xp 


Equating  Eqs.  (34)  and  (35),  solving  for  the  term  containing  integrals  for 


6yLy  and  substituting  into  Eq.  (33)  we  have 


xb  -  tb  tb  -  xb 

6j(6y)  =  /  (ctyt)6y|  dx  +  /  Xyxx6yx|  dt 
xQ  tQ  tQ  xQ 


tb  -  Xb  Xb  tb  Xb  - 

+  /  tUyx  -  Uyxx)x]5y|  +  &*yx6y|  }dt  +  /  /  SyQdtdx 

to  Xo  Xp  tQ  Xo 


tb  Xb 

-  /  /  {ayt6yt  +  Xyxx6yxx  +  £yx6yx  +  £*yx6yxH(x-xp) }dtdx  =0  (36) 

t-o  xo 


This  is  the  key  equation  which  uses  variational  principle  in  solving  a 


mixed  initial  and  boundary  value  problem  for  a  fourth  order  partial 


differential  equation. 


DISCUSSION  OF  THE  VARIATIONAL  EQUATION 


Let  us  discuss  the  various  terms  in  Eq.  (36),  the  variational  equation 


for  the  beam  problem,  into  three  pajrts  as  follows. 


(1)  The  initial  conditions  of  the  original  variables  are  given  and 


variations  of  the  adjoints  at  the  far  end  are  zero.  The  first  term  in  Eq. 


(36)  contains  the  product  of  y^Sy  evaluated  at  the  initial  condition  yto^yo 


and  at  the  final  condition  ytb^yb*  Since  the  value  of  yb  are  known  as  given 


by  Eqs.  (27)  and  (29),  Syb  =  0.  That  is,  the  variations  of  the  adjoint 


variable  at  the  far  end  are  zero. 
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(2)  The  boundary  conditions  of  the  original  variables  and  variation  of 
the  adjoints  can  be  determined.  The  second  through  fourth  terms  are  the 
boundary  terras  involving  the  variations  6y  and  6yx  and  the  variables  yx ,  yxx , 
and  yxxx  at  both  boundaries.  For  a  beam  the  end  conditions  can  be  expressed 
as 


Fixed  End 

'■C 

ii 

o 

y  =  0 

6y  =  0 

yx  =  o 

yx  =  0 

6yx  =  0 

Hinged  End 

"<J 

11 

O 

^  1 

11 

o 

<5y  =  0 

yxx  =  o 

yxx  =  o 

6yxx  =  1 

Guided  End 

O 

11 

X 

1 

X 

11 

o 

6yx  =  0 

yxxx  =  o 

yxxx  -  o 

6yxxx  = 

Free  End 

yxx  =  o 

yxx  -  o 

5yXx  ■  1 

yxxx  =  o 

yxxx  =  o 

6 yxxx  = 

The  variations  in  the  adjoint  variables  shown  in  the  last  column  coincide  to 
the  same  end  conditions  in  the  original  variables  given  in  the  first  column, 
whether  it  is  on  the  left  or  the  right  boundary.  It  is  noted  that  the  third 
partial  derivatives  can  be  evaluated  at  the  boundaries. 

(3)  Interior  region  -  The  last  two  terras  give  the  interior  where  the 
forcing  function  Q,  the  adjoint-variations  6y,  Syt ,  6yx,  and  6yxx  and  the 
variables  y^,  yx,  and  yxx  are  shown.  No  third  order  partial  of  y  with  respect 
to  x  is  present.  Thus  the  variables  that  are  needed  for  the  computation  are 
y»  yt,  yx»  yxt>  yxx>  an<i  yxxt*  This  requires  a  c2  continuity  in  the  spatial 
direction  and  a  c1  continuity  in  the  time  domain. 
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TRANSFORMATION  OF  COORDINATES 


The  integral  signs  in  Eq.  (36)  can  be  converted  into  summation  signs  if 


discrete  intervals  for  integration  are  used. 

We  may  take 

some 

scale  factors 

to  nondimensionalize  the  problem 

by  giving 

rr 

O 

II 

o 

cb 

=  1 

0  <  t  <  1 

(37) 

o 

II 

o 

Xb 

=  1 

0  <  x  <  1 

(38) 

Moreover,  Eq.  (36)  can  be  discretized 

by  letting 

5  =  Ht  -  i+1 

0  < 

k  <  i 

i  =  1,2,. 

..  ,H 

(39) 

n  =  Kx  -  j+1 

0  < 

n  <  1 

j  =  1,2,. 

..  ,K 

(40) 

where  H  and  K  are  number  of  intervals 

for  t  and  x  respectively. 

Thus  the 

partial  derivatives  are 

yt 

8y 

3t 

-  H  87  - 

n 

Hy  5 

(41) 

yx 

8y 

8x 

dy 

=  K  —  = 

Kyn 

(42) 

yxx  = 

82y 

Sx*  “ 

3yx 

:  K - - 

3n 

K2ynn 

(43) 

yxxx  = 

8  3y 
8x3  = 

aYxx 

K - = 

3n 

=  K3yrmn 

(44) 

Use  of  Eqs.  (36)  through  Eq.  (44)  then  leads  to 
0  =  Sj(Sy) 


K  i  _  tb  1 

-  i  \  )  ] 6yU,j )  |  -  dq 

3-1  0 


t0  K 


H  i  Xb  l 

+  I  f  UKyn  -  (XK3ynn) n] 6yC± , j ) |  -  d£ 

i=l  0 


xQ  H 
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+  l  l]  (^K2yTin)K6yn(i»j)|Xb  -  dC 

i=l  0 


x0  H 


+  l  /  A*Ky„6y(i»J)|  -  d? 

i=p  0 


xp  H 


K  1  H  i  _  ,  , 

+  I  I  (  l  f  6yC 1 , j ) Q  -  d£}-  dn 
j=l  0  i=1  0  H  K 


K  1  H  !  _  _  .  1 

'  \  ^0  {J,  L  faH2y5(i,j)6y5(i,j)  +  ^KVnnSypp  +  *K2yn6yn]  -  d£}  -  dn 
j=l  u  i=l  u  H  K 

K  1  H  i  _  l  l 

~  l  J  {  l  j  [^*K2yT16ynH(x-xp)]  -  d£}  -  dn  =  0  (45) 

i=p  0  i=1  0  H  H  K 


GRID  SYSTEMS 

The  (24x1)  vector  has  a  grid  of  four  (6x1)  vectors  through 

Y4(i»j),  thus 

Y(i,j)  =  {[Yi(i»j)]T[Y2(i,j)]TIY3(i,j)]T[Y4(i,j)]T  (46) 

Each  of  the  (6x1)  vectors  has  six  components  consisting  of  the  function,  its 
first  and  second  partials  in  spatial  directions,  and  its  mixed  partials  in 
space  and  time. 


Yl 


(i.j)  = 


y(Ci,nj) 

y(5i,ij+i) 

yc(5i.Dj) 

Y3^»j  )  = 

y^^i.nj+i) 

yn(5i,nj ) 

yn(5i»hj+i) 

ych^i.hj ) 

yCn(?i»hj+i) 

Ynn^i.hj ) 

Ynn^i  *hj+i) 

yCnn^i.Dj) 

y?nn^i*nj+l 
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Y2 


U,j)  = 


yUi+l.rij) 

yUl+l ,nj+i) 

y^i+l.nj) 

yn^i+l»nj> 

Y4U,j)  = 

y^Ui+i.nj+i) 

yn(5i+i»nj+i) 

y^n^i+l»T1j) 

y^n^i+i  »nj+i> 

ynn(^i+l  ,rij  ) 

ynn^i+i^j+i) 

y^nn^i+l  ,r|j  ) 

y^nn^i+i’^+i) 

(47) 

If  we  increase  the  row  index  from  i  to  i+1 ,  then  the  grid  point  shifts 


down  by  one  step  and  the  following  holds 

y^(i+l,j)  =  =  Y4(^»3)  (48) 

If  we  increase  the  column  index  from  j  to  j+1  then  the  grid  point  shifts  to 

the  right  by  one  step  and  one  obtains 

Yl(i,j+1)  =  Y3(i>j>  Y2(i>J+1)  =  (49) 

The  following  diagram  shows  the  relationship  of  the  grid  system. 


Yl(i»j) 

Y3(i,j)  =  YjCM+I) 

Y3(i»j  +  D 

[yU.j)] 

[Y(i,j+l)] 

Y4U»j+l) 

|  | 

yjd+l.j) 

i’(i+l,j)  =  y)u+1  »j+l) 

Y3U+I  ,j+l) 

[Y(i+1  >j>] 

[Y(i+l  ,j+l)] 

Y2(i+1  >j  ) 

Y^(i+l,j)  =  Y2O+I  >j+l) 

Y4(i+l,j+l) 
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SPLINE  FUNCTION 


We  may  express  the  variables  and  6y(i>J)  in  Eq.  (45)  in  terms  of 

the  (1x24)  spline  function  aT(£,n)  and  the  (24x1)  node  point  function 
as  follows. 

y(io>(£,n)  =  aT(^,n)Y(i*J)  (50) 

where 

aT(5,n)  =  {[a1(^,n)]T  [a2(£,n)]T  [a3(£,n)]T  [a4(S,n)]T  (51) 

and 

6y(i,j)(5,n)  =  aT(5,n)6Y(1»j)  (52) 

A  typical  term  for  a  product  can  be  written  as 

6yU>3)y(i»j)  =  [  6Y< 1  *.1 )  ]Ta(5,n)aTU,n)Y(l  >3  >  (53) 

CONCLUSION 

A  bilinear  form  of  the  original  and  adjoint  variable  is  employed  in 
determining  the  coefficients  of  the  variations  of  the  functions  and  their 
first  derivatives.  There  is  no  term  involving  the  variations  of  any  higher 
derivatives  than  second.  In  solving  mixed  boundary  and  initial  value  problems 
of  a  fourth  order  partial  differential  equation  using  spline  functions,  the 
computation  may  be  simplified  considerably  if  the  variable  in  time  can  be 
truncated  into  arbitrary  sections.  The  entire  problem  is  divided  into  several 
strips  of  distinct  time  intervals,  each  strip  containing  mostly  the  boundary 
value  problem. 

The  variational  principle  for  spatial  and  temporal  problems  with  boundary 
and  initial  conditions  have  been  investigated.  This  variational  principle  is 
very  general  in  scope  and  can  be  applied  to  many  linear  partial  differential 
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equations.  The  principle  is  applicable  if  the  bilinear  concomitant  is 
identically  zero.  This  leads  to  the  requirement  that  a  set  of  end  conditions 
for  the  adjoint  systems  must  be  found  to  satisfy  this  condition.  Otherwise 
the  variational  principle  as  stated  may  not  be  applicable. 

The  beam  equation  (with  one  dimensional  spatial  direction)  satisfy  these 
variational  principles.  For  future  work  the  analytic  solution  of  these 
equations  using  finite  element  method  will  be  studied.  The  assembly  of  the 
elements  of  the  matrices  involved  in  the  formulation  will  be  demonstrated. 

The  stability  problem  in  numerical  solutions  on  these  equations  will  also  be 

investigated.  This  lays  the  foundation  for  the  gun  dynamics  problem  to  be 
studied  in  the  future. 
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APPENDIX 


THE  VARIATIONAL  PRINCIPLE 

A  more  accurate  estimate  can  be  made  by  constructing  a  variational 
principle. 5  By  using  the  adjoint  variable  y  as  a  Lagrange  multiply  we  have 

J[y»y]  “  <Qy>  +  <y,(QhLy)> 

=  <Q,y>  +  <y,Q>  +  <y ,Ly>  (Al) 

In  order  that  J  be  a  variational  principle  the  following  requirements  must  be 
satisfied . 

(a)  J  is  stationary  about  the  function  ys  which  satisfies  the  following 
relation 

Lys  “  “Q  (A2) 

(b)  The  stationary  value  of  J  deduced  from  Eqs.  (2)  through  (5)  is 

J[y,y]  =  <Q,ys>  ♦  <Q,ya>  (A3) 

where  ya  is  the  actual  solution.  Consider  first  the  stationarity  of  J  by 
taking  the  variation  of  Eq.  (Al) 

6J  =  <Q,<5y>  +  <6y,Q>  +  <6y,Ly>  +  <y,L6y> 

=  <<$y,(Ly+Q)>  +  <<Sy,(Ly+Q)> 

-  <<5y,Ly>  +  <y,L6y>  (A4) 

We  will  make  an  effort  later  to  impose  certain  conditions  in  order  that 
the  following  equality  holds: 

<y,LSy>  =  <6y,Ly>  (A5) 

where  L  is  the  adjoint  operator. 


5Stacey,  W.  M.  Jr.,  Variational  Methods  in  Nuclear  Reactor  Physics,  Academic 
Press,  1974.  “  "  ' 
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By  combining  Eqs.  (A4)  and  (A5)  one  obtains 

6 J  =  <6y,(Ly+Q)>  +  <6y,(Ly+Q)>  -  0  (A6) 

Since  the  variations  6y  and  6y  are  arbitrary  it  leads  to  the  requirement  that 
the  stationary  values  ys  and  ys  must  satisfy 


Lys  =  -Q 

(A7) 

Lys  =  -Q 

(A8) 

Since  Eq.  (A7)  is  the  same  as  Eq.  (A2)  therefore,  J  is  stationary  about  the 
function  ys. 

Equation  (A8)  is  the  adjoint  equation  in  terms  of  the  adjoint  operator, 
L,  the  adjoint  variable  y,  and  the  adjoint  forcing  function  Q. 

It  is  noted  that  6J  in  Eq.  (A6)  vanishes  and  is  independent  of  .the 
arbitrary  variations  6y  and  6y.  By  using  6J  one  can  claim  that  the  estimate 
is  very  accurate  and  free  from  the  arbitrary  variations. 

Using  the  relationship  in  Eq.  (A7)  the  stationary  value  of  J  from  Eq. 
(Al)  is 

J[ys»ysl  =  <Q,ys>  +  <ys»Q>  +  <ys»Lys>  =  <Q»ys>  (A9) 

Since  J  is  stationary  and  6J  0,  then 

<Q,ys>  -  <Q,ya>  (aio) 

which  is  the  requirement  given  in  Eq.  (A3). 

It  is  noted  that  Eq.  (A6)  contains  no  boundary  terms  to  be  satisfied. 
This  bears  an  important  point  in  the  future  discussion. 
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